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Real scalar field models incorporating asymmetric double well potentials will decay to the state 
of lowest energy. While the eventual nature of the system can be discerned, the determination of 
the dynamics of the bubble wall provides many difficulties. In the present study we investigate 
numerically the evolution of spherically symmetric expanding bubbles coupled to a thermal bath in 
3 + 1 dimensions. A Markovian Langevin equation is employed to describe the interaction between 
bubble and bath. We find the shape and velocity of the wall to be independent of temperature, yet 
extremely sensitive to both asymmetry and viscosity. 
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I. INTRODUCTION 

The behavior of fields during weak first order phase 
transitions has been a source of great interest. Several 
studies have shown either a need for changes to the cur¬ 
rent theory of thermal phase transitions or the inclusion 
of carefully calculated corollaries to the present models 
10-H. Sub-critical bubbles and oscillon solutions may 
give rise to phase mixing in which these fluctuations in 
the smooth homogeneous background can have signifi¬ 
cant dynamical effects. Thus in the limit of weak first 
order phase transitions, the assumption of small ampli¬ 
tude fluctuations may need revision. For stronger transi¬ 
tions, the fluctuations are expected to diminish in impor¬ 
tance validating the use of standard homogeneous nucle- 
ation theory. Yet, in the realm of the expanding bubble, 
precious few analytic and numerical studies have been 
performed to determine the effect of the thermal bath on 
the expansion of the bubbles With this in mind, we 
numerically examine expanding bubbles in 3 + 1 dimen¬ 
sions using a Langevin equation in hopes of gaining either 
qualitative or quantitative results that describe their be¬ 
havior. 

Symmetry restoration at various stages in the Uni¬ 
verse’s life has become an important and powerful con¬ 
cept in modern cosmology. Several examples of the rele¬ 
vance of phase transitions are given by inflation 0, the 
electroweak phase transition Q, and the quark hadron 
phase transition 0. In particular, the phase transition 
that occured when the universe cooled to a critical tem¬ 


perature of about 300 GeV broke the symmetry of the 
weak and electromagnetic interactions. The breakdown 
of SU(2) X U(l) which may have created the baryon 
asymmetry we see today is expected to be of first order 
1^. This satisfies the third of Sakharov’s conditions for 
baryogenesis, namely that no thermal equilibrium exist 



At the electroweak scale, the rate of expansion of the 
universe is such that thermal equilibrium is maintained 
|pd| . To create an out of equilibrium process, a bubble 
of true vacuum appears within the false vacuum state. 
As the bubble expands, baryogenesis takes place in the 
neighborhood of the boundary of the bubble, the bubble 
wall. Thus, the detailed behavior of the bubble wall is of 
great importance in the electroweak phase transition. 

The nucleation of bubbles in the context of field the¬ 
ory has by now a long history. Coleman has shown that 
at zero temperature, the transition to the true vacuum 
takes place via the quantum nucleation of bubbles . 
This process corresponds to a generalization of barrier 
penetration in quantum mechanics. The probability for 
bubble nucleation at zero temperature is found by deter¬ 
mining the “bounce” configuration from the Euclidean 
equation of motion subject to the boundary conditions 
of restricting the field to the metastable minimum at 
Ie = —oo, the global minimum at Ie = 0, and the 
metastable minimum at = +cx). The bounce solu¬ 
tion is used to compute the Euclidean action, from which 
bubble nucleation probability per unit volume is deter¬ 
mined, r = A exp {—Se)- The dominant contribution to 
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the tunneling rate conies from the solution to the equa¬ 
tion of motion with the least action. These solutions are 
0(4) symmetric, satisfying a 0(4) symmetric Euclidean 
equation of motion. The nucleating bubbles were shown 
to quickly accelerate to the speed of light. 

Linde expanded Coleman’s ideas to account for finite- 
temperature effects in the nucleation of bubbles [ p^ . 
Temperature effects can be included by formal substi¬ 
tution of the Euclidean time to include temperature, im¬ 
position of periodic boundary conditions on the field, and 
integration in the “new” Euclidean time from 0 to 1 
For high temperatures, the kinetic term in the action 
grows large. Yet, as solutions which minimize the ac¬ 
tion are sought, the time independent configurations will 
dominate. The 0(3) Euclidean equation can be solved 
to find the finite temperature bounce. The rate of ther¬ 
mal nucleation of critical bubbles per unit volume was 
found to include an “activation energy” term represent¬ 
ing the free energy barrier the system must overcome for 
the transition to occur. The thermal fluctuations were 
shown to be a purely classical effect. Thus the previous 
results of the bounce solutions also apply to the finite- 
temperature case with different configurations for differ¬ 
ent temperatures. 

Recent numerical studies have shown that the actual 
first order phase transition from false to true vacuum in 
the presence of a heat bath may be a great deal more 
complicated than indicated from the picture of a finite 
temperature bounce solution. For weak first order phase 
transitions, phase mixing through sub-critical bubbles 
may effectively restore the symmetry Q . Long-lived sub- 
critical bubbles, oscillons, may act as nucleation sites for 
critical bubbles, thus speeding the transition rate [^,||. 
Although the importance of these effects are still some¬ 
what in dispute [|l5| , it is clear that significant deviations 
from the standard model of bubble nucleation may oc¬ 
cur and that carefully considered numerical calculations 
must be performed in order to determine the behavior of 
the system. 

Previous studies of thermal nucleation of bubbles have 
dealt with systems in lower dimensions. Work in 1-1-1 
dimensions has primarily focused on the nucleation of 
kink-antikink pairs II- The behavior that was exam¬ 
ined was the kink density and lifetime for a wide range 
of viscosities and temperatures. The numerical results 
showed good agreement with the theoretical expecta¬ 
tions. In 2-1-1 dimensions, studies have focused on the 
validity of theoretical calculation of the nucleation bar¬ 
rier when compared to numerical results |p^ and the ef¬ 
fect of fluctuations on the structure and velocity of the 
expanding bubble wall Q. 

In this study, we will investigate the behavior of a 3 -I-1 
dimensional bubble expanding in a thermal background. 
To insure uncontaminated results, the field is localized in 
the metastable vacuum with the heat bath acting until a 
bubble is formed. Although this requires longer compu¬ 


tational times, it provides results which occur naturally 
from the equations of motion without need for an arti¬ 
ficial contrivance. We hope to shed light on two vital 
characteristics of the transition from our study as the 
asymmetry of the potential, the bath’s viscosity, and the 
bath’s temperature are changed: 1. the shape of the 
bubble wall, and, 2. the speed at which the bubble wall 
travels. 

The rest of this work is organized as follows. In the 
next section, we briefly review the properties that give 
rise to expanding bubbles and the bounce solutions found 
by Coleman ||^. In Section III, we discuss the effects of 
coupling the system to a thermal bath by generalizing the 
equation of motion found in section II using a Langevin 
form. In Section IV, we present the numerical results. In 
Section V, we summarize our results, pointing to possible 
directions for future work. 


II. EXPANDING BUBBLES 


The behavior of a real scalar field follows from the form 
of the action. 


S = 



VW 


( 1 ) 


Through use of Hamilton’s principle, the classical equa¬ 
tion of motion for the field is found to be 


W 




dV[^) 

d(j) 


( 2 ) 


To determine the energy of the field, we integrate in space 
over the stress-energy tensor. 



l(V^)2 + u(^) 


( 3 ) 


where the dot refers to partial differentiation with respect 
to time. Thus, by specifying the potential, we will ob¬ 
tain an explicit equation of motion for the system and a 
determination of the energy of the field. 

In the early universe, and, particularly, in the case of 
the electroweak phase transition, the approximate ex¬ 
pression for the effective potential in the high temper¬ 
ature limit to one loop order is Q 

V(</., T) = D{T^ - - ET(I)^ + (4) 

where 
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b solve the equation of motion for our system, we 
|oduce three boundary conditions. The field initial 
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starts at rest, 


and 2/r^ = m|^, Inas = 21 n 47 r — 27 ~ 3.91, lna_F = 
2 In TT — 27 ~ 1.14. Analyzing the temperature depen¬ 
dence of equation we find that for large temperatures 
the potential exhibits a single unique minimum at 0 = 0 . 
As the system evolves and cools, a non-global second 
minimum forms at a temperature of 


v/1-9£;2/8AtiT> 

With further cooling, degenerate minima of the potential 
form at the critical temperature given by 


When T < Tc a, new global minimum is formed. At 
T = Tq, the potential’s global minimum is located at 

3ETo 

For T < Tc, first order phase transitions are possible 
for the system described by V{cj),T). The cubic term in 
the potential is responsible for barrier formation and the 
asymmetry between stable and metastable states. The 
strength of the asymmetry depends upon ET. To sim¬ 
plify the expression for equation ^ while retaining the 
parameters of physical interest and allow for possible ap¬ 
plicability to other similar first order scenarios, the po¬ 
tential is written as 


(5) 


4>(p,o) = 0 . 

To avoid an undefined value at the origin, the spatial 
derivative is fixed 


d$(0, t ) 

dp 


= 0 . 


Eventually the system will nucleate a bubble of true 
phase in the false one. Thus, at spatial infinity we expect 
the field to be in the false vacuum state, 

$(p ^ oo,t) = $0 ■ 



FIG. 1. Asymmetric double well potential with a = 2.6. 


The variable a will be a measure of the asymmetry be¬ 
tween the two minima. 

To cast the system into a more tractable form, two 
further steps can be taken to simplify the equation of 
motion. By imposing spherical symmetry, we effectively 
reduce this problem to one dimension. This allows for 
quicker computational manipulation and a clearer un¬ 
derstanding of what is driving the system’s behavior. By 
neglecting azimuthal and longitudinal components, this 
study will serve to provide bounds on the behavior of the 
expanding bubble. Introducing the dimensionless vari¬ 
ables p = rm, T = tm, $ = —d>, and a = the 
nonlinear equation of motion becomes 


Using the dimensionless variables, the characteristics 
of the potential can be determined. The minima of the 
potential are at 


$ = 0 


and 



For a 


, the potential is degenerate. 


As shown in 


figure 1 , when A > y f, ‘h-i- becomes the stable vacuum 
and the possibility for the creation of expanding bub¬ 
bles exists. The classical turning point of the potential 
is found at 
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III. BOILING BUBBLES 
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The zero temperature bounce and the finite tem¬ 
perature bounce solutions are needed in calculating 
the bubble nucleation rate. To calculate the bounce, we 
write a Euclidean version of equation ^ through use of 
t —ix^ and — |x|2 ^ —(x°)2 — jxp = —\xe\^- 

The boundary conditions restrict the field to be local¬ 
ized in one of the minima of the asymmetric double well 
potential at x° = ± 00 , hence the name “bounce”. The 
least action Euclidean solution will contribute most to 
the nucleation rate, having 0(4) symmetry. Recalling 
that a connection between the path integral formulation 
of quantum field theory and the statistical mechanics par¬ 
tition function can be found by the formal substitution 
of x° —> /3 with P = defining the functions on a 
domain of length (3 and periodically connected in time, 
the finite temperature action can be determined. The 
functional integral for the action now contains a kinetic 
term of quadratic order in temperature. At high tem¬ 
perature, this term can be neglected as we are looking 
for the solution which minimizes the action. In other 
words, the least action solution in finite temperature will 
be time-independent and have 0(3) symmetry. The 0(3) 
spherically symmetric equation of motion is 


(P<^ 2 _ 9U(4>) 

(ip2 p dp 5$ 


( 8 ) 


For the asymmetric double well potential, no closed- 
form solution exists to equation ||. However, approxi¬ 
mate solutions can be determined if one restricts the in¬ 
vestigation to asymmetries that are small compared to 
the barrier height of the potential. For the dimension¬ 
less form of the potential this amounts to a ^ In 
this limit, the volume energy dominates the surface en¬ 
ergy in equation ^ only when the radii of the bubbles 
are large. Additionally, the transition between the vacua 
in 4>(p) occurs only in a small interval, thus giving rise 
to the name “thin-wall” approximation . As the first 
order derivative term in equation ^ is only appreciable in 
a small transition region and p is large, we may approx¬ 
imate our equation of motion by 


02$ 


= (j) - 




(9) 


Solving this equation for $(p) we obtain the shape of the 
bubble in the thin-wall limit 


V2 


I — tanh 


p — R 


( 10 ) 


with R the bubble radius. This solution will be of use 
later in our investigation when we compare the numeri¬ 
cally generated bubble wall to the shape expected from 
a slightly asymmetric double well potential. 


To include the effects of a thermal background on the 
evolution of the field, a generalized Langevin equation 
can be formulated and solved numerically. Employing an 
analogy between classical Brownian motion and quantum 
field theory, we find an equation of motion by adding a 
viscosity and a noise term to our original equation of mo¬ 
tion, equation |[ Assuming a Markovian heat bath, the 
viscosity coefficient is related to the noise ^ (x, t) by the 
fluctuation-dissipation theorem 

(^(x,t)^(x',t')) = 2-fT6{t-t')S^{x-x') . (II) 


The coupling constants in this system are treated as free 
parameters. Although more complicated nonlocal forms 
of the Langevin equation could be used |Q , the increased 
complexity in both the equation of motion and the dy¬ 
namics do not warrant such an approach until the present 
situation is better understood. Indeed, the agreement 
between numerical studies and theoretical predictions of 
nucleation in I-l-I dimensions lends support to the use of 
an additive noise and a Markovian bath 

The system is restricted to spherical symmetry. As the 
thermal bath is provided with no symmetries, the proper 
evaluation of the system should contain all spatial di¬ 
mensions. Yet the advantages provided by evaluating 
the fully dimensional case is questionable in light of the 
adversity this would cause. The computer simulations of 
a non-symmetric system would be expected to have sig¬ 
nificantly longer run times. More importantly, the dom¬ 
inant perturbations due to the thermal fluctuations are 
expected to be along the radial direction. By allowing 
only variations in the radial direction, a clearly defined 
bubble wall velocity can be determined. For these rea¬ 
sons, the values obtained here should be considered as an 
approximate bound on the velocities of expanding (3-1-1)- 
dimensional bubbles. Writing the fluctuation-dissipation 
theorem in terms of spherical symmetry, we obtain 

(l(P: T)iip\ t')) = - t')6{p - p') , (12) 

where 9 = —T is the dimensionless temperature. The 
spherically symmetric dimensionless Langevin equation 
is written as 


92$ , 52$ 

i9r2 ^ dr dp^ 


2 9$ 9 - 

-k a$2 _ 

P op 

(13) 


where y = ^ and ^ are the dimensionless viscos¬ 

ity and noise respectively. 
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Equation ^ is numerically solved by a finite difference 
routine second order accurate in time and fourth order 
accurate in space. The step sizes employed in the sim¬ 
ulation were determined rather crudely by specifying a 
numerically obtained bounce solution as the initial field 
configuration in the 7 = 0 case and ensuring that this 
system behaved as expected for a sufficiently long time 
(r > 6000). Satisfactory results were found in these situ¬ 
ations if a spatial step of 5p = 10 “^ and a temporal step 
oi 6r = 2 X 10“^ were used. 

The validity of equations and can be tested by 
comparing numerical results and theoretical predictions 
for the system subject to a quadratic potential. Given 
sufficient time, the thermal evolution of the field in the 
potential will reach equilibrium. Since the system is de¬ 
scribed by a spatial one dimensional lattice of N grid 
points, we can measure the energy per degree of freedom 
E/N. From basic statistical mechanics, we know that 
once the system is in equilibrium the equipartition theo¬ 
rem will give E/N = 6/2. The numerical results for this 
situation do show that the equipartition theorem holds 
for our system over a wide range of parameters [See Fig¬ 
ure 2 of Ref. HI]. 

While we have explicitly formulated a system whose 
dynamics are given by the Langevin equation, many hy- 
drodynamical studies have shed light upon the phase 
transition ||20|-p^. Whereas in the Langevin descrip¬ 
tion the particles scattering off the bubble wall provide 
a damping proportional to the field’s change in time, a 
hydrodynamic description treats the wall as a combus¬ 
tion front. The results of the hydrodynamic work have 
shown various characteristics of expansion which are not 
manifest using the Langevin equation. Instabilities have 
been found in the background fluid and the bubble wall 
which may dominate the dynamics 1^ . The velocity 
of the bubble wall has also been found to be dependent 
upon the thermal conductivity of the background fluid 
and upon temperature inhomogeneities arising from the 
release of latent heat These considerations, while 
interesting and important, are not investigated in this 
work. The dynamics of the bubble we obtain here come 
from the simplified description of the field subject to a 
potential interacting with the heat bath by equation 

IV. NUMERICAL RESULTS 

Equation ^ behaves in a manner consistent with ex¬ 
pectations in the basic case of equilibrium thermodynam¬ 
ics. Computational grid step sizes have been found that 
reproduce behavior present in specific aspects of nucle- 
ation. Our investigation into the effects of a thermal bath 
on the shape and velocity of an expanding bubble wall 
may now begin. 

For the field to be trapped in the metastable state and 
eventually decay by bubble nucleation to the global vac¬ 


uum, <5 > y|. As shown in figure 2 , the bubble is cre¬ 
ated from the background field fluctuations. Eventually 
the fluctuations are large enough to form a critical bub¬ 
ble, a bubble that is large enough to complete the phase 
transition and grow. To distinguish bubble formation 
from spurious background fluctuations computationally, 
we examine the field searching for values greater than the 
turning point. Once a segment of length 10 in p has been 
found whose members satisfy this criterion, it is marked 
as an expanding bubble. 

The shape of the bubble wall in figure 2 is reminis¬ 
cent of the hyperbolic tangent thin-wall configurations. 
Generalizing equation |^, we obtain 


= ^[l-tanh(p/)] (14) 


where is the global minimum and f~^ specifies the 
thickness of the wall. In figure 3a and 3b, equation 
is used to fit expanding walls generated in baths with 
a = 2.8, 7 = 1.0, T = 3.0 and a = 3.1, 7 = 2.0, T = 1.0 
respectively. While the smoothness of the wall is altered 
when temperature is varied as asymmetry and viscosity 
remain constant, the bubble wall’s mean shape does not 
change. 



FIG. 2. Snapshots of bubble evolution from the back¬ 
ground field. The potential has an asymmetry a = 3.1 and 
the heat bath is characterized by 7 = 1.0 and T = 1.0. The 
interval between snapshots is Ar = 10“^. 
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Figure 3a 
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Figure 3b 
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FIG. 3. Snapshots of expanding bubbles with equation 
providing best ht given by equation 14 for (a) a = 2.8, 
7 = 1.0, T = 3.0 and (b) a = 3.1, 7 = 2.0, T = 1.0. 

The dependence of the wall thickness on the asymme¬ 
try of the potential and viscosity of the thermal bath can 
be found numerically by holding one quantity constant 
while varying the other. For 7 = 1.0 and a ranging be¬ 
tween 2.4 and 3.4, hgure 4 shows that the data for / can 
be described by a linear equation 

f(a)=Aa-B (15) 

with A = 1.21 and B = 2.20. In figure 5, the numerical 
results of / as asymmetry is held constant and viscosity 
varies is shown to be well fit by the equation 

f{j)=Ar^+B (16) 


with A = 0.512 and B = 0.969. 



2.4 2.6 2.8 3 3.2 3.4 


FIG. 4. /(d) vs. asymmetry for constant viscosity 

7 = 1.0. The equation providing best fit is Td — B where 
A = 1.21 and B = 2.20. 
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FIG. 5. /(t) vs. viscosity for constant asymmetry 

a = 3.1. The equation providing best fit is -|- B where 

A = 0.512 and B = 0.969. 


To ascertain the kinematics of the bubble wall, a com¬ 
putational technique for determining the size of the bub¬ 
ble must be developed. As the turning point interpolates 
between the two vacua, the size of the bubble can be 
defined as that radial distance which corresponds to the 
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field configuration value being equal to the turning point. 
Recording the bubble size at a particular time allows one 
to find the wall velocity through the equation Vhw = ^ 
where Ap is the difference between the bubble size at 
successive times and Ar is the difference between the 
corresponding times. In order to limit the already in¬ 
tensive use of computational resources employed in this 
study, the simulations were terminated once the bubble 
reached a size of 80. Several trials with longer cutoff sizes 
were performed and found to yield no additional benefit 
for the parameters probed. 

The development of the bubble wall velocity through 
time is shown in figure 6. The graph clearly indicates a 
region in which the bubble wall is accelerating, eventually 
reaching a terminal velocity. Since the region of acceler¬ 
ation is only of a limited duration, we will ignore this, 
deferring a detailed analysis of the feature to later stud¬ 
ies. Additionally, the restriction imposed of identifying 
an expanding bubble by a region of true vacuum phase 
of size 10 (in p) and greater in most cases obscures the 
accelerating wall behavior. Thus, although the data ob¬ 
tained exhibits an acceleration in the start of a bubble’s 
life, our computational methods restrict us to a purely 
qualitative analysis of this phenomena. 
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ble. For added certainty, a second method of smoothing 
the velocity as a function of time data using a Savitzky- 
Golay type filtering scheme was employed concurrently 
. By examining the graph, the velocity could be visu¬ 
ally discerned. Both methods produced values in agree¬ 
ment with each other. 

As in the case of the wall shape, the velocity of the wall 
was found to be invariant to the temperatures of the bath 
probed here (T = 10“'’’ to T = 10). To insure that this is 
not an artifact of some statistical dependence in the noise 
term, all computational simulations were given different 
random number seeds and a reliable random number gen¬ 
erator used PH . The generator was also checked for any 
correlations and none were found. The parameters affect¬ 
ing the velocity were the bath viscosity and the potential 
asymmetry. 

The bubble wall terminal velocity’s dependence on 
asymmetry of the potential for constant viscosity 7 = 1.0 
is shown in figure 7. The best fit equation to the numer¬ 
ical data is provided by 


’^hw (c^) 


1 

1 -I- Aa + Bct^ 


+ 1 


(17) 


with A = 1.61 and B = —1.16. 



FIG. 6 . Bubble wall velocity vs. time for a = 3.1, 7 = 0.5, 
and T — 0.5. The dotted line indicates the constant bubble 
wall velocity 0.935. 


FIG. 7. Bubble wall velocity vs. asymmetry for constant 
viscosity (7 = 1.0) and temperature (T = 1.0). The equa¬ 
tion providing best fit is + 1 where A = 1.61 and 

B = -1.16. 


The value for the velocity of the bubble wall is deter¬ 
mined numerically through use of a “windowing” least 
squares method. The data over which a least squares fit 
for a linear equation was performed was adjustable so as 
to find those portions which had a slope < 10“®. 

Once this criterion had been satisfied, the intercept was 
interpreted as the terminal expansion velocity of the bub- 


In figure 8, the asymmetry was fixed to a = 3.1 while 
the viscosity was allowed to vary. The velocity’s depen¬ 
dence is well approximated by the expression 

vUl)= +-P ( 18 ) 

with A = 0.949, B = 0.275, C = 1.62, and D = 0.0584. 
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FIG. 8 . Bubble wall velocity vs. viscosity for constant 
asymmetry (a = 3.1) and temperature (T = 1.0). The 
equation providing best fit is + D where A = 0.949, 

B = 0.275, C = 1.62, and D = 0.0&4. 


try which tends to entice the system back to the expected 
results. For velocity, the dependence on asymmetry is 
approximated by equation 0 - As the asymmetry of the 
potential increases in the presence of a heat bath, the ve¬ 
locity will approach the expected value of unity, the speed 
of light. Yet, as is shown by equation 18, the viscosity 
slows the bubble wall. Here the equation behaves as ex¬ 
pected in the limit 7 —> 0 giving light speed velocities for 
the expanding bubble wall. If equation 18 is carried to 
large viscosity, the unsettling behaviour of velocities for 
all viscosities exists. In fact, letting j 00 one finds 
a suggestive lower bound on expansion speeds of 0.0584. 
This would indicate that viscosity alone can never pre¬ 
vent the bubble from growing. Yet, taking this result 
too seriously may lead to error due to the lack of numer¬ 
ical trials of sufficiently large viscosities performed and 
the methodology employed. It must be borne in mind 
that our investigation has focused on only growing bub¬ 
bles which evolve naturally from the false vacuum. For 
large viscosities, the possibility of finding well-defined, 
expanding bubbles becomes questionable. This possibil¬ 
ity in turn may lead to a discontinuity in an extended 
investigation of velocity vs. time, eventually leading to 
the expected result of zero velocity. 


V. CONCLUSIONS 

We investigated the shape and velocity of expanding 
bubble walls in the presence of a thermal bath. The 
dynamics of the system were modeled by a generalized 
Langevin equation with a Markovian bath and an asym¬ 
metric double well potential. The bubbles were gener¬ 
ated naturally from the background thermal field with¬ 
out need for specifying an initial configuration. In the 
absence of the thermal bath, expanding bubbles eventu¬ 
ally reach a terminal velocity of the speed of light and, 
in the thin wall limit, will have a profile corresponding 
to a hyperbolic tangent shape. 

The presence of a thermal background strongly affects 
the shape and velocity of the wall. While the relevant 
aspects of the system were invariant under changing tem¬ 
peratures, asymmetry and viscosity altered the dynamics 
of the expanding bubbles. The shape of the wall is well 
approximated by the thin wall shape with a thickness pa¬ 
rameter that is dependent on both viscosity and asym¬ 
metry. The inverse thickness of the bubble for changing 
asymmetries is given by equation In the presence of 
a heat bath, as the asymmetry increases, the thickness 
of the bubble wall decreases. For changing viscosity, the 
shape is given by equation |^. The thickness of the bub¬ 
ble wall is found to increase as the viscosity increases. 
Unlike the expected results of increasing wall thickness 
as asymmetry grows for no thermal background [ pd| , the 
results here indicate that a relationship exists between 
viscosity which thickens bubbles and increasing asymme- 


Our investigation focused on spherically symmetric so¬ 
lutions to a Langevin equation with a Markovian bath. 
Yet, there remain several questions whose answers would 
provide much new information on the subject of first or¬ 
der phase transitions. To model truly real-world phe¬ 
nomena, ideally one would employ a fully dimensional 
equation of motion. The prospects of doing this present 
serious challenges not only in the amount of present-day 
computational resources needed but also in such concepts 
as clearly defining a bubble wall velocity. An interesting 
extension of this study appears to be in the determi¬ 
nation of the acceleration of the walls from bubble cre¬ 
ation to the time of terminal velocity seen in the first 
few moments of figure 6 . More fundamental to the in¬ 
vestigation of interacting field theories, it is possible that 
the thermal background will furnish more complicated 
properties such as non-additive couplings and/or nonlo¬ 
cal correlations in time and space [p^ . The consequences 
of these nonlocal effects are still questionable, but it is 
clear that their inclusion into the dynamics of the system 
by coupling to more general thermal baths may have very 
significant effects on the shape and speed of the bubble 
walls. Of more fundamental theoretical importance is the 
possible connection between the electroweak phase tran¬ 
sition and the evolution of the system according to the 
Langevin equation. While the methods employed here 
are strongly suggestive of the electroweak phase transi¬ 
tion, the work of proving the connection between the two 
is left as an open question. 
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